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Abstract. We investigate the ground states of classical Heisenberg spin systems 
which have point group symmetry. Examples are the regular polygons (spin rings) 
and the seven quasi-regular polyhedra including the five Platonic solids. For these 
examples, ground states with special properties, e. g. coplanarity or symmetry, can be 
completely enumerated using group-theoretical methods. For systems having coplanar 
(anti-) ground states with vanishing total spin we also calculate the smallest and 
largest energies of all states having a given total spin S. We find that these extremal 
energies depend quadratically on S and prove that, under certain assumptions, this 
happens only for systems with coplanar S — ground states. For general systems the 
corresponding parabolas represent lower and upper bounds for the energy values. This 
provides strong support and clarifies the conditions for the so-called rotational band 
structure hypothesis which has been numerically established for many quantum spin 
systems. 
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1. Introduction 

The study of small interacting spin systems is not only of theoretical interest but also of 
importance for the experimental investigation of recently synthetized magnetic molecules 
0,0,0,0- An exact calculation of the thermal expectation values of the relevant quan- 
tum observables and other quantities such as correlation functions is possible only in 
very few cases. Given this situation, a classical treatment yields a first approximation 
for individual spins s with s ^> 1 which is often astonishingly good, as well as bounds for 
the exact quantum values (c. f. Berezin/Lieb-inequalities). It is also valuable as a guide 
for developing approximation schemes for attacking the quantum theoretical problem. 

One of the fundamental characteristics of a system is its ground state(s) and the 
corresponding ground state energy. The problem of determining the exact classical 
ground state has been considered for a long time, see e. g. @,||, but has only been 
solved for a few special systems. These include arrays of interacting spins occupying the 
sites of a regular iV-polygon (spin rings) with nearest-neighbor coupling constant J > 
(antiferromagnetic in our convention). The corresponding ground state energy is 

E min = -2JNs\ (even N), (1) 

corresponding to an alternating configuration of spin orientations in the iV-polygon. 
The analogous result 

E min = 2JNs 2 cos((N- l)n/N), (odd N) (2) 

does not appear to be in the literature, although it is easily established, see (fT3|). An- 
other class of systems where the classical ground states are easily determined are the 
so-called bi-partite systems, which can be divided into two subsets A and B such that 
there are only interactions (with a positive coupling constant J) between spins a e A 
and b G B. In this case any anti-parallel configuration between A and B will minimize 
the energy. The regular polygons with even N considered above and the cube are simple 
examples of bi-partite systems. 

Recently, the exact classical ground state of a system consisting of iV = 30 spins 
occupying the vertices of an icosidodecahedron [7j],|| has been determined [[|. The 
method employed there can be applied to any system which can be decomposed into 
a set of triangles without common edges and which is 3-colorable. Examples are the 
octahedron (as mentioned in ||) and the cuboctahedron (this is obtained by joining 
the mid-points of the edges of a cube with their nearest neighbours, see figure |6|). All 
these systems have coplanar ground states with S = which can be modified to remain 
ground states in the presence of a magnetic field 7i, the corresponding energy depend- 
ing quadratically on TC, c. f. |§. This yields a definite prediction for the magnetisation 
versus magnetic field plot at low temperatures: It should have a constant slope until a 
critical magnetic field Ti c is applied, and for H > H c the magnetisation is saturated. 
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In the present paper we will derive a further consequence for systems having 
coplanar ground states with 5 = 0. The same formal equation as that characterizing the 
ground state in the presence of a magnetic field occurs for the ground state subject to the 
extra constraint that the square of the total spin assumes a given value. Generalizing the 
solution in |9| we obtain a whole family of states assuming the minimal energy E m i n (S). 
It turns out that for these systems the function E min (S) is quadratic in 5 of the form 

E m US) = J -^j^S 2 +j mm N. (3) 

Here j is the row sum (J = 'Y1 IU J^ v ) of the adjacency matrix J defined below. J is 
chosen in such a way that j does not depend on /x, and j m m is the smallest eigenvalue 
of J. Conversely, the existence of a lower parabola (§) implies the coplanarity of the 
ground state, see theorem || for details. Analogous results hold for the maximal eigen- 
value E max (S) of systems with 5 = coplanar anti-ground states. Hence these systems 
show an exact rotational band structure (RBS) which has been conjectured to apply 
for a general class of quantum spin systems jl(J. Thus our findings strongly support 
that conjecture in the sense that the considered class of spin systems has exact RBS- 
parabolas in the classical limit s — > oo. On the other hand we have found classes of 
systems which deviate from the exact RBS-parabolas, especially for small values of 5, 
see section 6. Hence our work may be viewed as a first step towards an understanding of 
the conditions for the occurence of RBS-parabolas. We will call systems, which exactly 
attain their lower RBS-parabola fl3|), "parabolic". For general systems the correspond- 
ing parabolas represent lower and upper bounds for the energy values, as will be proved 
in section 4. Sometimes these bounding parabolas will be very good approximations of 
the boundaries of the energy spectrum even for non-parabolic systems. 

Other questions arise in connection with spin systems having coplanar 5 = ground 
states. Has every classical spin system coplanar ground states or even coplanar 5 = 
ground states? If not, what is the criterion for having coplanar ground states? Are there 
systems with both coplanar and non-coplanar ground states? We do not have complete 
answers to these questions, but we can show the following: 

• In general the spin triangle with arbitrary coupling constants J\,J2,Jz has only 
coplanar ground states with 5 > 0. 

• A spin system which we call "pentagonal star" has a non-coplanar ground state 
with 5 > 0. The existence of coplanar states with the same energy can be excluded 
on the basis of numerical calculations. 

• For certain systems the existence of coplanar, symmetric ground states can be 
excluded. By "symmetric" we mean roughly that the state has the same symmetry 
as the spin system itself (for details see section 5). 

• Other systems, such as the dodecahedron and the icosahedron, have symmetric 
non-coplanar ground states and apparently no coplanar ones. These ground states 
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can be geometrically vizualized in terms of what is known as the "great stellated 
dodecahedron" and the "great icosahedron" , respectively. 

• The tetrahedron has a symmetric non-coplanar ground state with S = and a 
variety of other non-symmetric ground states including coplanar ones, all of the 
same energy. 

• The same is true for other systems with full permutational symmetry ("iV- 
pantahedron"), except that these systems do not have symmetric ground states 
for N > 4. 

• The cuboctahedron has a lower RBS-parabola but interestingly it has two 
symmetric ground states with S = 0, a coplanar one and a non-coplanar one. 
The latter corresponds to a degenerate stellated figure which we have not found in 
the literature and which is composed of four Stars of David, see figure |[ 

Thus we are faced with a surprising variety of possibilities of ground states of differ- 
ent types. The main open problem is to rigorously prove the existence of spin systems 
without coplanar ground states. It appears that quantitative statements can be made 
mainly if symmetric states are concerned, since in this case we may use standard meth- 
ods of group theory, in particular representation theory of the point groups. 

The article is organized as follows: Section 2 contains the general definitions and 
results which are independent of symmetry assumptions. The main result is that systems 
with coplanar S = ground states have RBS-parabolas (theorem 1) and its converse 
(theorem 2). These parabolas represent lower and upper energy bounds even for non- 
parabolic systems, as proven in section 4. In section 3 we address questions of extension 
of ground states to larger systems and "higher order frustration" . Section 5 introduces 
symmetry assumptions and the definition and simple properties of symmetric states. 
This machinery is applied in section 6 to the investigation of particular examples. We 
consider the general triangle which is generically non-parabolic as well as classes of 
symmetric spin systems including the quasi-regular polyhedra. We construct symmetric 
ground states from certain irreducible representations of the respective point groups. 
The article closes with a summary. 

2. Notations, definitions, and general results 

We consider as the phase space of the classical spin system the iV-fold product of the 
unit sphere S 2 . Instead of using canonical coordinates we will denote a state by a se- 
quence s of unit vectors s M , n — 1, . . . , N with components s l , i — 1, 2, 3. The total spin 
is S = J2n s v with components S\i = 1,2,3. If not indicated otherwise, the abstract 
letters . . . will always be (upper) indices and there is no danger of confusion with 
the square of a vector, e. g. S 2 . We will use bracketed indices to denote by, say, a?W the 
vector with the components x l , similarly for matrices. For example, the n x m matrix 
consists of n row vectors and of m column vectors A^j. 
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If not mentioned otherwise, the Hamilton function (or "energy" ) will be of the form 
H (s) = ^ J ^ s v ' s »- ( 4 ) 

If a magnetic field 7i is to be included we add a Zeeman term and obtain 

H h = H -h-S, (5) 



where h = ys(s + l)gnsH contains the common combination of the spectroscopic 
splitting factor g and the Bohr magneton and is scaled with the function a/ s(s + 1) 
of the corresponding spin quantum number s. 

Note that the exchange parameters J^ v are not uniquely determined by the 
Hamiltonian H via Different choices of the leading to the same H will 

be referred to as different "gauges". We will adopt the following gauges. First, the 
antisymmetric part of J does not occur in the Hamiltonian @. Hence we will follow 
common practice and choose J^ u = J Vfi . Thus the J^ v can be considered as the entries 
of a real symmetric matrix J. Second, since s M • s M = 1 we may choose arbitrary diagonal 
elements J MM without changing H , as long as their sum vanishes, TrJ = 0. The usual 
gauge chosen throughout the literature is J MM = 0, fj, — 1, . . . , N, which will be called 
the "zero gauge". In this article, however, we will choose another gauge, to be called 
the "homogeneous gauge" , which is defined by the condition that the row sums 

will be independent of \i. Note that the eigenvalues of J may non-trivially depend on 
the gauge. We found that the eigenvalues of J are only relevant for energy estimates if 
the row sum of J is constant. This would restrict the applicability of large parts of our 
theory if we adopt the zero gauge. However, by introducing the homogeneous gauge we 
can apply our results to a very general class of systems. 

The quantity 

Nj EE £ J„u (7) 

fJbV 

is gauge-independent. If exchange parameters satisfying J^ IU = J Vii are given in the 
zero gauge, the corresponding parameters J^ u in the homogeneous gauge are obtained 
as follows: 

■V = Jfu, for /Ji^v, (8) 
Jy,n — 3 -Jfj.- (9) 

It follows that 
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A spin system is called antiferromagnetic (AF) iff all J^ u > for // ^ v. A system 
is called connected iff its spin sites cannot be decomposed into two disjoint classes, say 
A and B, such that J^ v = if ji G A and v G B. 

In the special case where all exchange parameters J lxVl fi 7^ v are or J, the system 
can be essentially represented by its graph T. A graph T = (V, £) consists of a set of 
"vertices" V and a set of "edges" 8 C [V] 2 . Here [V] 2 denotes the set of subsets of V with 
exactly 2 elements. In our case, the vertices are the spin sites, V = {1, . . . , iV}, and the 
edges represent interacting pairs of spins. For spin systems it is appropriate to consider 
only graphs without loops and multiple edges, as we do in this article, following, e. g. , 

There are different ways to graphically represent a spin state s M . For specific spin 
systems, e. g. magnetic molecules, the spin sites // G V are embedded into the physical 
3-space and represented by vectors r^. One way would be to attach the vectors as 
small arrows to the sites given by r /t . This is done in figure |7| below. However, when 
using this representation it is difficult to visualize the structure of the in spin space. 
Hence we will mostly employ another method of representation: We draw the N vertices 
in spin space according to the unit vectors and join them with lines according to the 
edges of the original spin graph T. Thus we obtain a graph isomorphic to T which 
contains the additional information about the considered spin configuration. 



In graph theory the set of edges is often represented by a matrix called the 
adjacency matrix. We will use this name also for the matrix J in the general case 
of different exchange parameters. There exists an extended literature about the con- 



nection between the structure of a graph and the spectrum of J, see, for example fl2 
where also applications in chemistry and physics are mentioned. 

The graph T is called complete if £ = [V] 2 . We will also call the corresponding spin 
system where any two spins interact with equal strength, a pantahedron. 



Being symmetrical, J has a complete set of N real, ordered eigenvalues j m in, • • • , jmax 
which are counted according to their multiplicity. One of them is the row sum j with 
1 = 77^(1) 1, • • • , 1) as the corresponding eigenvector. For connected AF-systems, j will 
be non-degenerate and equals the largest eigenvalue, j = j max > 0, by the theorem of 
Frobenius-Perron. It will be convenient to rearrange the indices such that j = jo and 
the remaining eigenvalues ji, ■ ■ ■ ,jN-i ar e ordered according to their size. Note that 
for AF systems, ji = j min since j min = j > contradicts Tr J = 0. However, since it 
may happen that j max = j we will write j ma xi for the largest eigenvalue of J other than 
j, counted once, hence Jn-i = jmaxi for AF systems. Sums over a = 1, . . . , N — 1 but 
excluding a = will be denoted by ^ . We will denote the a-th normalized eigenvector 
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J2 J ^ C «=J* C ^ ^« = 0,...,iV-l, (11) 

and 

^C?C£ = 5 a/3 , Qt,/3 = 0,...,JV-1, (12) 

where we also allow for the possibility to choose complex eigenvectors. 

We note in passing that J in the homogeneous gauge is essentially the Hamilton 
operator of the corresponding quantum spin system, restricted to the subspace of total 
magnetic quantum number M = Ns — 1. Thus there is an unexpected connection be- 
tween "weakly symmetric" classical states to be defined below and quantum eigenstates 
of the Hamiltonian lying in the mentioned subspace. 

The equations of motion resulting from the Hamiltonian (|J) are 

d 



dt 




S K = I 7 , J ta> 8 v I X V ( 13 ) 



Using = J Vfl one can immediately show that the total spin vector S is a constant of 
motion. 

We now will define various kinds of special states. A stationary state will be one 
with is/j, = 0, for all /i = 1, . . . , N. According to (0), this is equivalent to 

V 

for some real numbers k m , \x = 1, . . . , N. Equation (|HD can also be viewed as the solution 
of the problem to seek states with vanishing variation of the quadratic form 

^ ^ J^lU S fl ' s u (^5) 

subject to the constraints s M • s M = 1, \i = 1, . . . , N. The k m then appear as Lagrange 
parameters of this variational problem. 



Geometrically, (14) means that for a stationary state each spin vector is propor- 
tional to the "center of mass" of its neighbors. As for general Hamiltonian systems, 
the stationary states are just the critical points of the Hamilton function, i. e. those 
points where the gradient of Hq vanishes. In particular, the states s with minimal 
H (s) = E min or maximal H (s) = E max are stationary. The former will be called 
ground states, and the latter anti-ground states. For AF systems the anti-ground state 
is always of the form = e for all // = 1, . . . , N, i. e. all spins are aligned parallel to an 
arbitrary unit vector e (ferromagnetic ordering). In this case j max = j and E max = Nj. 
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Moreover, we will consider the relative ground states, (relative anti-ground states), 
which are defined as the states of minimal (maximal) energy among all states satisfying 
(^2 u s u ) 2 = S 2 . Let E min (S) denote the energies of the relative ground states, and 
Emax(S) that of the relative anti-ground states. Again, the relative ground states and 
anti-ground states are among the solutions of the variational problem corresponding to 
(P^D with the additional constraint s u ) 2 = S 2 > 0. If the extra Lagrange parameter 
is called %■> we obtain the condition, analogous to (|14|), 

^2 J ^ s " = K ^ s v + X& for S > 0. (16) 

V 

The equations of motion fllTf ) then imply that for states satisfying 

^-s^ = XS x s M , (17) 

i. e. , for 5 > these states are no longer stationary but are precessing around the total 
spin vector. This is completely analogous to the precession of stationary states of H in 
the presence of a magnetic field. 



16) we have 



The condition 5 > in fllq ) is necessary in order to apply the method of Lagrange 
parameters. 5 = has to be excluded since it would correspond to a critical point of 
the constraining function g(s) = s u ) 2 — S 2 , see for example JT3|. To cover also the 
case 5 = we consider three constraining equations written in vector form 

$>„ = So, (18) 

V 

where Sq is a fixed vector. Then the above-mentioned problem does not occur. We have 
three Lagrange parameters which will be written as the components of a vector L. The 
resulting equations are 

JjwSy = + L. (19) 

v 

Taking the vector product of (|19|) xs M and summing over fi yields o = S x L. If 5^0 
this is equivalent to L = %5 and we recover (|16[). For 5 = 0, however, (|T6|) has to be 



replaced by (19) 



States satisfying (|T9"D will be called weakly stationary states, since they are station- 
ary in a rotating frame. Hence relative (anti-) ground states are weakly stationary. 

There exists a simple mechanical model for stationary states: Consider a system 
of iV rigid, massless rods of unit length fixed at the same point ("pendula"). These 
rods may be represented by unit vectors s^. Between two rods s M and s u we consider 
springs satisfying Hooke's law with spring constants 2J flu . Further there is a constant 
("gravitational") force F = Fe. By the cosine theorem, — s u \ 2 = 2(1 — s M • s u ), 
hence the potential energy of our mechanical model will be 

V(s) = )^2 J ^\ S » - s »\ 2 ~ F ■ S = N 3 J ^ s » -Sv-F-S- (20) 
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The equilibrium points of this mechanical system will satisfy the zero-force equation 



where 2/t^s^ are essentially the constraining forces exerted on the rigid rods. Moreover, 
for equilibrium S will be proportional to the constant force, say F = —2xS, since oth- 
erwise there would be a first-order variation of the potential energy when performing 
infinitesimal uniform rotations of the system. It follows that the equations ([H]) and ( [L6[) 
are formally identical. Summarizing, there is a 1 : 1 correspondence between the weakly 
stationary spin states s and the equilibrium states s of the considered mechanical sys- 
tems, if the value of the force F runs through all real numbers. This mechanical model 
may help the reader to vizualize some of the weakly stationary states to be considered 
in this article. 

We can easily obtain the following bounds for the energy of stationary states s: 

Lemma 1 Nj min < H Q (s) = ^k,, < E max < Nj max . Moreover, for AF 

systems and S = we have H (s) < Nj maxi . 

Proof: H (s) = J^s v ■ s M = k m s m • = ^ V Further we have J2^u J ^ s ^ 4 
s n < 3maxYli^ s iJ. ' s n = Nj maxi and analogously for the lower bound. Since (anti-) 
ground states are stationary, these bounds also hold for E min and E max . The upper 
energy bound for S = and AF systems follows from J^s v ■ « M < Jmaxi E M s v ' s m> 
since for S = o the vectors s^, % = 1,2,3 are orthogonal to the eigenvector l corre- 
sponding to the eigenvalue j. ■ 

The respective bounds of lemma 1 are attained if there exist stationary states where 
all = j min or k m = j max , or = j maxi in the AF case for S = 0. At present, we 
cannot show the existence of such states under general conditions but this can be done 
for some specific examples. We will call any stationary state with the property that 
the values of all k m of fll4l) are independent of /x, a weakly symmetric state. If we view 
the set of components of a state s as an iV x 3-matrix, the definition of "weakly 
symmetric" could be rephrased in the way that all rows of s are eigenvectors of J for the 
same eigenvalue j. Recall that for stationary states each spin vector is proportional 
to the "center of mass" of the neighbouring spins. For weakly symmetric states the con- 
stant of proportionality will be the same for all spins s M , which motivates the wording. 
Obviously, weakly symmetric (anti-)ground states belong to the eigenvalue j m i n , (jmax)- 
If one takes any three eigenvectors of J with the same eigenvalue to form the 3 rows of 
a matrix, one will not get automatically a weakly symmetric state, since the iV columns 
need not be unit vectors. In section 5 we will present a construction procedure which, 
under certain assumptions, yields (weakly) symmetric states. 




(21) 



Since the energy Hq(s) is invariant w. r. t. rotations/reflections in spin space the 
above-defined classes of weakly symmetric states, (weakly) stationary states, or of (rel- 
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ative) ground states are also invariant under rotations/reflections. 

Another simple property of weakly symmetric ground states is the following: 

Lemma 2 Every weakly symmetric state has vanishing total spin if the corresponding 
eigenvalue 7 is different from the row sum j . In particular, every weakly symmetric 
ground state of an AF system has vanishing total spin. 

Proof: We obtain J ^ s » = 3 s v = 3 s and 12^ J ^ s » = s n = I s - Hence 
J = j or S = o. If s is a weakly symmetric ground state of an AF system with j — jmim 
we have j min < j as above. ■ 



We will call a state s collinear if all its spin vectors are (anti-)parallel to a given 
vector e, i. e. = ±e. Due to equation (14) , every collinear state is stationary. Further 
we will call a state s coplanar if all its spin vectors ly in a plane i. e. , if there exists 
a non-zero 3- vector n such that n ■ s M = for all \i = 1, . . . , N. In section 6 we will 
show that systems such as the regular polygon, the pantahedron and the quasi-regular 
cuboctahedron and icosidodecahedron admit coplanar weakly symmetric 5* = ground 
states. These ground states are important since one can construct from them a whole 
family of relative ground states, which show a rotational band structure. 

Theorem 1 (i) Let s be a coplanar weakly symmetric state with S = corresponding 
to an eigenvalue J of J and perpendicular to a unit vector n. Then the following 
family of states s(S), < S < N has total spin length S: 



C2 a 

s^S) = ^l-—s, + -n, p=l,...,N. (22) 

(ii) If, moreover, s is a ground state then s(S) will be a relative ground state for all 
< S < N. In this case, 

E mm (S) = J ~^ mm S 2 + 3mm N. (23) 

(Hi) The analogous case for s being a coplanar weakly symmetric relative anti-ground 
state for S = 0: Then s(S) will be a relative anti-ground state and 

E max (S) = tlJpi s 2 + j maxi N. (24) 

Proof: 

(i) Obviously, s(S) is a state with total spin length S. 

(ii) We obtain J2„ = \/l - jpjminS„ + §jn. This is of the form (0) if we 
choose x = ■ > ~ 3 ™ in . In order to prove that s is a relative ground state, we note 
that s is a weakly symmetric ground state w. r. t. the modified adjacency matrix 
J^ v = J^ u — x, since j min is also the smallest eigenvalue of J. The calculation of 
E m in(S) is straightforward. 

(iii) This case is largely analogous, but j min has to be replaced by j ma xi, since s has 
S = and cannot be the total anti-ground state. 
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Systems which satisfy equation will be called parabolic systems. Hence the 



essential claim of theorem [1] is that systems with a coplanar S = ground state will be 
parabolic. 

We will now prove a theorem which can be viewed as the converse of theorem 1, 
in so far it shows that coplanar ground states necessarily appear for parabolic systems, 
if certain technical assumptions are satisfied. We expect that these assumptions hold 
under fairly general conditions but will not dwell upon this question further. 

Theorem 2 Consider a connected AF system and let t i— > s(t) be a smooth curve for 
< t < e of relative ground states such that the limits t — > of s(t) and its derivatives 
up to second order exist, in particular ^f(O) ^ for S 2 (t) = s M (i)) 2 . Moreover, 
let s(0) be a weakly symmetric ground state with S = and assume parabolicity, 
E m i n (S) = 3 ~ 3 ™ in S 2 + Nj min , at least for the interval covered by < t < e. 
Then s(0) will be coplanar. 

Proof: We will indicate the limit t — > of a function of t by skipping the argument, 
e. g. s M = 8^(0). Differentiation with respect to t will be indicated by an overdot. Since 
s M (t) ■ Sfj,(t) — 1, differentiation yields 

s M ■ = 0, s M ■ + ■ = 0. (25) 

We may assume that S(t) = S(t)e where e is a constant unit vector. Being a weakly 
stationary relative ground state, s(t) satisfies fliTf ) for t > 0: 

V V 

with the limit value S = s v = o according to the assumptions of the theorem. But 
we cannot a priori assume that has a finite limit value for t — > 0. Thus we solve 



(ID for L{t) = S(t)x(t) and obtain 



( 5 ~ l-(e-^(t)) 2 ■ (27) 

If s is collinear, the proof is done. If s is not collinear, we find some \i such that 
(e ■ s^(t)) 2 ^ 1 for all < t < e. For this value of fi, (p7|) defines a smooth function of 
t, independent of /i, with a finite limit value L = L(0). Moreover, we can solve (p6| ) for 
K^(t), — 1, . . . , N and obtain 

= J ^ s ^(0 ' s mW " L W e ' (28) 

This shows that also ft M (t) is a smooth function of £ and has a finite limit for t — *■ 0. Since, 
by assumption, s is weakly symmetric this limit must be = j min and hence L = 0. 
Therefore the smooth function x(t) = for t > has a finite limit x = h m t-^oo j^j, 
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employing l'Hospital's rule and that S(0) ^ by assumption. 



Differentiating ( p6|) twice and taking the limit t — > yields 



and 



^ ] JfiuSf — '^a s a + + jmin 'S^ + 2^S + X<S\ 



Multipying (p6|) , (p9|) , and ( |30| ) with s M and summing over [i yields 

H = Jgu S fi ' s u = ^ K fi 4~ = iV j m i n , 

~ UV fl V 

(J.V 

jmin ^ ^ 4~ 4" ^ ^ ^ 4" jmin ^ ^ ^ ' s ^ 

= 4-x^ 2 , 

A 1 

using ( pop and (^). On the other hand, 

H(t) = J2^(t) + x(t)S\t), 

H{t) = M*) + XWS 2 (*) + 2 X (t)5(t)5(t), 

and 

// • -A.v4 

Comparison with (|35| ) yields ^ fy» = an d 
= 2 X S 2 . 

Hence, using the assumption of parabolicity and l'Hospital's rule, 



j jmin dH 



H 2 X S 2 



d 2 Q2 OC2 



N d(S 2 
Summing (^) over /i yields 

jS = ^2 k ^ s ^ + ^ min ^ + N X&> 



X- 
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whence, by ( (4Tf ) and the assumption S ^ 0, 

$>m«m = °- ( 42 ) 



We rewrite (|29|) in the form 



5J("V ~~ 3miJ>iu> ~ NxE fMU )s u = k^s^, (43) 

where E is the projector onto J's eigenvector l with constant entries = -h. We 
expand both sides of fl43| ) into real eigenvectors of J and obtain 

a 

s v =J2'-^-K^ + ±S. (45) 
Now ■ «m = 0, • ^ k„s„ = 0, and ^ C£C% = 5 a p imply 

=EE'^ • E'- 2 -^ = £' ( 46 ) 

ft a jj J@ J a Ja J 

Since j > ja for a > and connected AF systems, this is only possible if all K a = o. 
Hence, by ( flop s u = j^S and • iS 1 = 0. Thus s is coplanar. ■ 
An analogous theorem concerning E max (S) and coplanar anti-ground states can be 
proved similarly 

3. Extension of ground states 

We consider spin systems with J^ u = or J and which can thus essentially be 
characterized by their graph T = (V, £). A sub-graph (representing a sub-system) 



T = (V, £) is defined by the conditions V C V, £ C £ and £ C 



V 



Let 



H = 2J s ^ s - ( 47 ) 

be the Hamiltonian of the subsystem. Let s be any state of T and s be its restriction 
to T. If s happens to be a ground state of T, s need not be a ground state of T, but, of 
course, H (s) > E min . 

The inverse problem is the extension problem: Given a ground state s of T, can it 
be extended to a ground state s of T (such that s is the restriction of s)? In general, this 
will not be possible. For example, if T is the graph of a dimer with ground state s =TI, 
the extension of s to a global ground state is only possible for bi-partite systems. We 
have no general theory describing the obstacles against extension (which could be called 
frustration of higher order), but we can describe the typical situation where extension 



is possible. For an analogous treatment in the quantum case, see JT4 . 
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The situation we have in mind is one where the graph of the spin system Ti is 
decomposed into smaller graphs which are copies of some graph T 2 in such a way that 
each edge lies in exactly one of the copies. For example, an icosidodecahedron can be 
decomposed into 20 triangles with disjoint edges. We have local states by considering 
the ground states of the copies of T 2 . But these local states must fit together in order 
make it possible to construct a global ground state of IV For the example of decompo- 
sition into triangles the condition of fitting-together of local ground states is equivalent 
to the 3-colorability of the large system. The icosidodecahedron happens to satisfy this 
condition and hence possesses a ground state of the energy which is 20 times the ground 
state energy of the triangle. For general T 2 the colorability condition has to be replaced 
by the existence of a map from Fi to T 2 which maps edges onto edges. Such a map may 
be called a "graph homomorphism" . 

A graph homomorphism h : Fi — > F 2 is a pair of maps h = {h^,hs) such that 
hv '■ Vi — > V 2 , h.£ : Si — > S 2 and {/i, u} G £\ implies h £ {fi, u} = {/iy(/i), hv(v)}. The 
definition of a graph isomomorphism is analogous. 

Proposition 1 Consider two spin systems with graphs Fi and T 2 and the following 
properties 

(i) there exists a graph homomorphism h : F\ — > T 2 , 

(ii) there exists a decomposition into sub-graphs T\ = U«=i ^V, which is a disjoint 
decomposition with respect to edges, 

(Hi) the restriction of h to F^, : — > T 2 is a graph isomorphism for each 
/i — 1 . . . k. 

Moreover, let : V 2 — ► S 2 be a ground state ofT 2 . 
Then = h v o will be a ground state ofTi. 

Proof: According to the disjoint decomposition of Si, Hi will be a sum of terms 
which are each minimized by the state Hence s^ 1 ) is a ground state of IY ■ 

The above-mentioned construction of ground states for the icosidodecahedron, the 
octahedron and the cuboctahedron follows the description given in this proposition. 
Other examples show that the two given conditions of the proposition are essential: The 
dodecahedron is 3-colorable but it cannot be decomposed into edge-disjoint triangles. 
Actually its ground state energy is lower than that of the state obtained by the 3- 
coloring. A simpler example is the bi-partite 3-chain which can be mapped by a graph 
homomorphism onto the triangle but is not isomorphic to the triangle. On the other 
hand it is possible to connect 4 triangles such that their edges are still disjoint but 
form a tetrahedron which cannot be 3-colored. This is then a "higher-order frustrated" 
system with a ground state energy which is larger than 4 times the ground state energy 
of a triangle. 
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4. Bounding parabolas 

In this section we will show that the parabolas which occur in theorem |I| are general 
bounds for the extremal values E min (S) and E max (S), even if the system is not parabolic. 
The proof is completely analogous to the quantum case, see P5I. We recall that C ( a u) 
denotes the a-th normalized eigenvector of J and consider a transformation of the spin 
vectors analogous to the transformation onto the eigenbasis of J. Define 

Definition 1 T a = Y.^ s ^ and Q a = T a ■ T a a = 0, . . . , N - 1. 

The inverse transformation then yields 

«a. = Z! C '« T «' M=l,---,iV. (48) 



Especially, T = S/yN since a = corresponds to the eigenvector l. The following 
lemma follows directly from the definitions: 

Lemma 3 N = E,M 2 = E a Q a = j;S 2 + £' Q Q a . 

Our main result of this section is formulated in the following theorem: 

Theorem 3 The following inequality holds: 

J Jmin ry2 , ■ i\t ^ tt ^ J Jmaxi /-»2 »r 

^ S +JminN <H < — S +j maxi N . (49) 

Proof: We rewrite the Hamilton function in the following form and conclude 

H =J2 J ^cfr a ■ Tp = Y^^Qp = jf s2 + y E,'ipQp ( 5 °) 



^ jyd ' Jmin / j 



ft 



3 -.s 2 + 3mm (n- Is 2 ) = L^s 2 + 3mm N 



N J \ N J N 
using (|4"8"D, the positivity of Qp, and lemma ^|. The other inequality follows analogously! 



We note that the proof does not depend on the dimension of spin space, but the 
optimal energy bounds may depend on this dimension. The regular pentagon is an 
example of a system where the above bounding parabolas are assumed by 3-dimensional 
states, but not by 2-dimensional, i. e. , coplanar states, see section 6.3.1. 



5. Symmetric states 

In this section we will assume suitable symmetry properties of the spin system and cor- 
respondingly consider the notion of "symmetric states". We make use of some simple 
concepts and results of group theory, which may be found in many textbooks, e. g. |16| ]. 

Let J^y = or J and consider the graph T = (V, £) characterizing the spin system. 
Further let Q denote the symmetry group of the graph T. Hence Q consists of all 
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permutations of V which map edges onto edges. Q has a "natural" representation by 
real N x iV-matrices, by permuting the standard basis of M. N in the same way as the 
spin sites /i G V. We will denote the iV x iV-matrix representing any symmetry g G Q 
by g and by G the set of all N x iV-matrices obtained in this way. Obviously, 

Lemma 4 The adjacency matrix J commutes with all g G G. 

Hence 

Lemma 5 The eigenspaces of J split into orthogonal direct sums of irreducible 
subrepresentations of the natural representation of Q . 

The irreducible representations of the relevant groups are well known, see e. g. [ 16|| , 
and can easily be associated to the different (subspaces of) eigenspaces of J. 



For specific spin systems, e. g. magnetic molecules, the spin sites /i G V are 
embedded into the physical 3-space and the group Q could also be identified with one 
of the point groups, i. e. finite subgroups of 0(3, M). There are only a finite number of 
possibilities, see e. g. [fL6| , 3.1.2. and 3.1.3. The most interesting cases are 

• the dihedral group V n , 

• the tetrahedral group T, 

• the octahedral group 0, 

• the icosahedral group y, 

as well as the improper point groups attached to them. 



In what follows, we make the following general assumption: 

Assumption 1 All exchange parameters J^ v are or J, and the group Q operates 
transitively on V. 

The latter means that for each pair fi, v G V there is a permutation g G Q such that 
g{jj) = v. The two groups Q and 0(3, M) operate independently on the set of states: If 
two matrices g G G and R G 0(3, M) are given, we define their action on a state s by 
matrix multiplication with s which is again considered as an iV x 3-matrix: 

Definition 2 (R, g) • s = i?sg _1 . 

The transformations R G 0(3, M) will also be referred to as "rotations/reflections 
in spin space" . 

We have the following obvious results: 
Lemma 6 Let g G G and R G 0(3, M.) and s be a state. Then 

(i) s and i?sg _1 have the same total spin length S, 

(ii) If s is (weakly) stationary, then Rsg' 1 is (weakly) stationary, 
(Hi) If s is weakly symmetric, then i?sg _1 is weakly symmetric, 
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(iv) if s is a (anti-) ground state, then Rsg 1 is a (anti-) ground state. 

Having denned the action of the product group (9(3, R) x G on states, it is natural 
to consider the corresponding subgroup leaving a given state invariant: 

Definition 3 For any state s let Q s = {(R, g) G (9(3, R) x G|i?sg _1 = s}. Q s is called 
the symmetry group of s. 

Definition 4 A state s is called symmetric if the projection onto the second factor 
7r 2 : Q s — ► G is surjective, i. e. if for each g G G there exists a rotation/reflection 
R G 0(3, R) such that R^sg = s. 

One may thus say that for symmetric states s every permutation g G G of the 
vectors s M can be compensated by a suitable rotation/reflection in spin space. 

For example, if (V, £) is the graph of a regular iV-polygon, the coplanar symmetric 
states are in 1 : 1 -correspondence to the roots of unity, z^ = 1, namely 

= z£ = exp(i/zn27r/iV), n, /i — 1, . . . , N. (51) 

Here Q = V N , the dihedral group of order 27V. 

We have the following 

Lemma 7 Let be a non-vanishing N x 3-matrix such that for all g G G there exists 
some R G (9(3, R) such that sg = Rs. Then, for some suitable A G R, As will be a 
symmetric state. 

Proof: It remains to show that Xsfi will be unit vectors for all [i — 1, . . . , N. Fix some 
indices /i, v. Using assumption 1, we choose a g G G such that (sg)^ — s v — Rs^. Since 
R G (9(3, R) the two vectors s u and have the same length. Because fj,, v were arbi- 
trarily chosen, all vectors s M have the same length, say A -1 , which completes the proof. ■ 

If a state s is not collinear, then Rs = s implies R — 1 In this case for each 
g G 7r 2 (C? s ) there exists a unique R G (9(3, R) such that i? _1 sg = s. We will write 
R = p s (g). It is easily shown that the map g h- > p s (g) is a linear representation of the 
subgroup 7T2{Gs) C G. If s is collinear we may uniquely fix R by choosing R G {1, —1} 
and thus obtain a representation p s also in this case. 

For the case of symmetric s we have 7i 2 (Q s ) = G and hence p s will be an n- 
dimensional representation of G with n < 3. 

In the case of a symmetric ground state it is thus necessary that the eigenspace of 
J corresponding to the smallest eigenvalue j m i n contains an irreducible representation 
of G of dimension n < 3. If it does not, we can exclude symmetric ground states. 

For practical purposes of constructing (anti-)ground states it would be desirable to 
invert the process and to reconstruct s from a given n-dimensional subrepresentation 
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p of the natural representation of Q, n < 3 . We will now describe this procedure and 
consider, for sake of simplicity, only the case n — 3. Let S be an 3- dimensional subspace 
of ~R N left invariant by all g G G and sl\, i = 1,2,3 a basis of S. It follows that the 
three vectors 5^07*0)3^, i = 1, 2, 3 are linear combinations of the basis vectors, hence 

E^4 = E^)^> (52) 

or, in matrix notation, p(g)s = sg. Being a real representation of a finite group, p is 
equivalent to an orthogonal representation, i. e. p(g) = Tp(g)T~ l is orthogonal with 
some invertible 3 x 3-matrix T. It follows that 

p(g)(Ts) = (Ts)g, (53) 

which means that Ts is a symmetric state with px s = P- Note that the orthogonality 
of the representation is crucial in order to obtain a row (Ts)f p = 1,...,N of unit 
vectors, compare lemma Hence we have proven the following 

Proposition 2 Let p be a 3- dimensional subrepresentation of the natural representation 
of Q . Then there exists a 3- dimensional symmetrical state s such that p s = p. 

Further there holds 
Proposition 3 Each stationary symmetric state s is weakly symmetric. 
Proof: We will use a matrix notation and write k = diag(K±, . . . , Kn). Then the con- 



dition (|16D for stationarity may be written as sj = sk. For g 6 G we have k = 
g _1 Kg = diag(K e (i), . . . , « g (jv))- We conclude s/tg = sjg = sgj = p s (g)sj = p s (g)sK = 
sg(g _1 /tg) = p s (g)sK. Hence sk = sk, which is in components k^s^ = K g ( M )S M whence 
k^ is independent of p (using assumption 1). ■ 

The converse of proposition [3] is not true: There are ground states of the 
icosidodecahedron which are weakly symmetric, but not symmetric, see section 6.5. 
The condition of stationarity in proposition ||] can be replaced by another condition: 

Definition 5 Let Q s denote the subgroup of Q leaving p e V fixed and G(p) the matrix 
group generated by its natural representation. A symmetric state s will be called isotropic 
if p s (G(p)) contains at least one rotation with an angle a ^ 0,tt. 



Lemma 8 If s is isotropic and p s = p g > then s M = ±s' 



A 4 ' 



Proof: s M is invariant under all rotations/reflections p s (g),g G G(p). Since some of 
these is a rotation with an angle a ^ 0, it, the axis of rotation will be unique and p s = p s > 
implies s M = ±s' M . ■ 



Proposition 4 Each isotropic state is stationary, hence weakly symmetric. 
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Proof: Using the above matrix notation we obtain sjg = sgj = p s (g)sj hence sj is a 
symmetric (not normalized) state with p s j = p s . Using that s is isotropic and lemma 
(§) we conclude (sJJ) M = ±As M , hence s is stationary, and, by proposition ([|), weakly 
symmetric. ■ 



6. Examples 

We will mainly consider AF systems and set J > throughout this section if not 
mentioned otherwise. 

6.1. The general spin triangle 

We consider spin systems with N = 3 and general coupling coefficients Ji, J 2 and J 3 . 
The special symmetric case J\ = J 2 = is atypical and its properties are discussed 
later. Another special case is the 3-chain with J\ = 0, J 2 = J3, which is probably the 
simplest example of a non-parabolic system. Its investigation will be left for the reader. 

Let s u , v = 1, 2, 3, be any spin configuration and consider the 3 x 3-matrix S with 
coefficients 

Sfo, = s^-s u . (54) 

It is positive and its diagonal consists of l's. Conversely, any 3 x 3-matrix with these 
properties is of the form (0). This follows from the spectral theorem. These matrices 
can be written in the form 

^1 w v ^ 

S= w 1 u =[u,v,w], (55) 

L v u 1 l 

where u, v, w are real numbers subject to the constraint 

det S = 1 - {u 2 + v 2 + w 2 ) + 2uvw > 0, and w 2 , v 2 , w 2 < 1. (56) 

Hence the set V of the matrices of the form ( |54|) can be considered as a compact convex 
subset of 1R 3 . It has the form of an "inflated tetrahedron", see figure |l|. Since the map 

tt:V — ► 1R 2 (57) 
ir[u, v, w]= (3 + 2{u + v + w), J\U + J 2 v + J^w) (58) 

= (5 2 (s),^ (s)) (59) 

is affine, the set V 2 = ^[P} of all possible values of the energy vs. square of the 
total spin will be a compact convex set too. To simplify the geometry we will assume 
J\ + J 2 + J3 = 0. Then P 2 can be considered as the orthogonal projection of V onto a 
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suitable plane. If we introduce new coordinates in the u, v, w-space by 

u + v + w (60) 
JiU + J 2 v + J 3 w (61) 
(J 2 ~ Js)u + ( J 3 - Ji)v + (J x - J 2 )w (62) 

the orthogonal projection n is essentially the projection onto the a, e-plane. Hence the 
boundary of Vi can be obtained by solving the two equations 



a 



T 



detS 



0. 



VdetS 



0. 



(63) 
(64) 



/ 



\ h - -h 

The first equation ( |63| ) defines the boundary of V, the second one (|4]) is satisfied by 
those points of the boundary of V which are orthogonally projected onto the boundary 
of V%. By using computer algebra software it is straightforward to express (|63l) and ( p4T) 
as equations for a, e, r and to eliminate r between (]63f) and (164]). Yet the result for the 
general case is too complex to be reproduced here. We choose the special case 



Ji = 0, J 2 = 1, J 3 = -i 
and obtain the following equation for the boundary of V 2 

= S 2 {9 - 10S 2 + S 4 ) 2 + (-27 + 288S 2 + 18S 4 
- 8(3 + S 2 ) 2 e 4 + 16e 6 , 



24S 6 + S 8 )e 2 



(65) 



(66) 



where we have re-substituted for a. The actual boundary of V2 is only a part of 
the family of curves defined by (|66D and is shown in figure |2|. Here we have an exam- 
ple where the curves E min (S) and E max (S) can be explicitely calculated and are not 
parabolas, although they are close to their bounding parabolas which are represented 
as broken lines in figure 0. This is not in contradiction with theorem [1| since the ground 
state of this system is coplanar but has S > 0. 



We expect that similar figures also appear for the general system of iV spins. 
Parabolic systems are probably rare cases, some of which will be considered in the 
remaining subsections. 



6.2. Sn -symmetric systems 

Recall that the iV-pantahedron was defined as a system of N spins where each pair 
(fi, v) interacts with equal strength J > 0. These systems serve mainly as exactly solv- 
able model examples for the sake of illustration, but the dimer, the triangle, and the 
tetrahedron are also of experimental relevance pT7[],[]18||,[P-9y . The symmetry group of the 
iV-pantahedron is obviously the group Sn of all permutations of iV spin sites. 
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By completing squares the Hamiltonian can be written as 

H = J(S 2 -N). (67) 

Since the energy depends only on S, each state is a relative (anti-)ground state with 

E min (S) = E max (S) = J(S 2 - N). (68) 

For S = there exist symmetrical ground states only if iV < 4. These are the obvi- 
ous dimer, triangle, and tetrahedron configurations in spin space. The dimension dim of 
the irreducible representations of Sn is given by the number of standard Young tableaus 
(see [0) and hence dim = 1 or dim > N — 1. The one-dimensional representations are 
either the trivial one corresponding to the total anti-ground state or the one attaching 
to each permutation its sign, which is not generated by a symmetric state. Thus there 
are no symmetrical ground states for N > 4. For example, a cyclic permutation of 3 
spins leaving 2 other spins invariant cannot be compensated by a reflection/rotation in 
spin space. 

The coplanar iV-polygon configuration in spin space does however define a weakly 
symmetric ground state with = —1. Hence the construction of theorem 1 yields 
symmetrical relative ground states s(S). We will check our equation (^) for this case: 
Let E denote the N x iV-matrix with E^ v = 1/N for all /i, v. It is the matrix of the 
projector onto the one-dimensional subspace spanned by the vector l = -^(1, 1, • • • , 1)- 
The adjacency matrix for the iV-pantahedron is J = J(NK —1). Hence its eigenvalues 
are j = J(N — 1) with eigenvector l and jmm — — J with the (N — l)-dimensional 
eigenspace of vectors orthogonal to l. Thus j min = j ma xi an d the two functions E min (S) 
and E max (S) coincide in accordance with what has been said above. 



6.3. T>n- symmetric systems 



6.3.1. The N -polygon The iV-polygon has the dihedral group TV as symmetry group. 
We have already mentioned the set of S — 0, stationary, coplanar, symmetric states s 
given by 

a„ = z£, /i = 0,...,iV-l, (69) 

where z n = exp(in2n / N) , n — 1, . . . , N, is a root of unity, z% = 1. Indeed, these states 
are the complex eigenvectors of the corresponding adjacency matrix 



J 



/o 

1 




1 



1 









1 











1 





1 







1 



1 \ 




1 



(70) 



The eigenvalues are obtained by 

(j«)„ = j^+^r 1 ) 



2Jcos(27m/iV)s A 



(71) 
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and the corresponding energies are 

E n = 2JNcos(2im/N), n = 0, . . . , N - 1. (72) 

It follows that 

p N - j E N/2 = -2JN for N even, 

mm Jmm \ E {N . m = 2JNcos(n{N-l)/N) foriVodd 1 } 

and 

£ m a* = E = Nj = 2NJ. (74) 

This follows since the bounds of lemma 1 are attained by these states. For N = 5 
the ground state of the pentagon may be vizualized as a pentagram in spin space and 
analogously for other odd N. 

Moreover, 

2NJcos{2tt/N). (75) 
The relative ground states of theorem 1 have energies 

E (<!)-l 2J ( 2 W~ N ) foriVeven, 
mm[ } ~\ 2J((l-cos(vr^i))f + jVcos(7r^!)) foriVodd, 1 ' 

and 

E max (S) = 2J((1 - cos(27r/iV)^ + iVcos(2vr/iV)). (77) 
Both parabolas meet in the anti-ground state with E max = 2JN for S = N. 



6.3.2. Coplanar states of the pentagon The regular pentagon as a special case dealt 
with in the preceding subsection assumes its bounding parabolas. However, the extremal 
energies E min (S) and E max (S) can only be realized by non-coplanar states, except for 
S = and S = 5. This has not been rigorously proven but shown by numerical evidence, 
see figure [| The spectrum E versus S 2 realized by coplanar states is a subset of the 
full spectrum with concave boundaries E^fl anar (S) and E^ l x anar (S) . 

The permutation of spin sites (0 <-> 0, 1^2^4^3^1) leaves the 
Hamiltonian invariant and maps relative ground states onto anti-ground states, hence 
E c ^l anar (S) = S 2 - 5 - EZi anar {S). The (anti-) ground state with S = 1 and E = -7 
(resp. E = 3) is especially remarkable: Its state vectors s M occupy the vertices of 
a regular hexagon leaving one vertex free. The anti-ground state (E^ l ^ nar (l) = 3) 
satisfies equation (|T6D with \ = 0, hence it is a stationary state, not only weakly 
stationary as other relative anti-ground states. Actually, E%^ nar has a local minimum 
at S = 1, see figure [| For axisymmetric weakly stationary states in the neighborhood of 
the hexagonal state, equation (|T6| ) can be solved analytically by using computer algebra 
software and hence E and 5* 2 can be expressed as functions of a common parameter, 
although of forbidding complexity. The resulting curve turns out to be a part of the 
upper boundary of the coplanar spectrum and is displayed in figure ||, as well as the 
corresponding lower boundary part. 
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6.3.3. The pentagonal star We obtain the pentagonal star (N = 6) by joining the 5 
vertices of the pentagon with its mid-point, see figure f|. This system has T> 5 symmetry, 
but it is not XVsymmetric in our sense, since does not operate transitively on the 6 
spin sites. Nevertheless, we will discuss this example since it has the interesting prop- 
erty that its ground states are not coplanar and E min (S), E max (S) are only piecewise 
parabolic. 

If we add to the 10 edges of the pentagonal star the 5 edges of the corresponding 
pentagram we obtain the 15 edges of the 6-pantahedron. Hence its Hamiltonian reads 

H = S 2 — Q — H 5} (78) 

where is the Hamiltonian of the pentagram, which is the same as that for 
the pentagon, up to a suitable permutation of the spin sites. This shows that 
the configurations minimizing H for given S are exactly those which maximize H^, 
analogously for E max (S). The maximal values for H 5 are given by the parabola, compare 

E%l(S 5 ) = 2(1 ~ C ° 5 S(27r/5)) g 5 2 + 10cos( 27 r/5), (79) 

where S5 is the length of the total spin of the 5 vertices of the pentagon. Here S$ has to 
be chosen maximal for given S, that means that S$ = S + 1 for < S < 4 and S5 = 5 
for 4 < S < 6. This yields 

E mm {S) = I s2 _ 16 if 4 < 5 < 6 (80) 

Hence the total ground state is attained not for S = 0, but for S = 3 ~ 2 V ^ = 0.381966 . . . 
with E min = — 5 — 2\^E = —9.47214 . . .. The corresponding spin configuration is shown 
in figure |5]. 

It remains to show that no other, coplanar ground state exists. Employing the 
numerical result of the preceding subsection for the pentagon, E^ l x anar (S) < E max (S) 
for < S < 5, it is easily shown that all coplanar states with S = have an 

energy larger than E m i n . The minimal energy for all coplanar states appears to be 
Emfn anar = — 9, realized by a hexagonal ground state with 5 = 0, but we do not yet 
have a rigorous proof for this claim. 

Analogously one can show that 

^f(S + 4)(S - 1 + Vh) if 1<5<6 , . 



Emax^yS^j 



6.4. O- symmetric systems 



S 2 -16 + ^±±^ if < S < 1 



6.4.1. The cube As already mentioned the cube allows a bi-partition and hence 
possesses a ground state of the form 

8„ = (-l)"e, // = !,..., 8, (82) 
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if the vertices are suitably labelled. From this we obtain relative ground states with 

E mm (S) = J(^ 2 - 24), (83) 

since j = 3J,j m i n = —3 J. Apart from j and j m i n the other eigenvalues of J are ±J 
with three-fold degeneracy, respectively. The corresponding eigenspaces carry two in- 
equivalent irreducible representations of O, called F 1: F 2 , see [pjfl . According to the 
considerations in section 5 we may conjecture that these two representations are gener- 
ated by stationary, symmetric, non-coplanar states. 

Indeed, this can be directly verified for the following states: Let r^/j, = 1, ... ,8 
denote the unit vectors pointing to the vertices of the cube. The the states 

sW = tv, /i = l,...,8, (84) 

and 



(-i) = (_l)/y /!=!,..., 8, (85) 



have the desired properties. They are easily vizualized: is just the "cube in spin 
space" and s^" 1 ) the tetrahedron where each of the four distinct spin vectors is attached 
to pairs of vertices of the cube connected by space diagonals, is the anti-ground 
state for S = with energy E max (0) = 8J. Since it is not coplanar, theorem 1 is not 
directly applicable. However, it is possible to project the cube s^ 1 ) onto a square and, 
in a second step, to enlarge the square to a square of unit vectors. Since these are linear 
transformations, the resulting co-planar state s^' is weakly symmetric and has the same 
energy as before, namely 8J = E max (0). Now theorem 1 yields E max (S) = J(\S 2 + 8). 

6.4-2. The octahedron As noted in ||, the octahedron can be decomposed into four 
triangles with disjoint edges and it is 3-colorable. Hence it has coplanar ground states 
with S = which are obtained by extensions of the local ground states of the triangles. 
Since j = 4J and j m i n = —2.1 we again obtain the RBS-parabola 

E mm (S) = AS 2 - 12). (86) 

The eigenspace of J corresponding to j m i n is two-dimensional and carries a real, 
irreducible representation of O. This corresponds to the one-dimensional complex 
eigenspace of J spanned by the vector 

Sf, = z n ^), /i = l,..., 6, (87) 

where the 

z n = exp(m27r/3), n = 0, 1, 2, (88) 

are the 3rd roots of unity and fi i— > n{jj) denotes any 3-coloring of the octahedron. It 
follows that the state s is symmetrical. 
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Table 1. Eigenvalues j n of the adjacency matrix JJ of the cuboctahedron (1st column) 
together with their degeneracy (2nd column). In the 3rd column the irreducible 
representations of O are indicated which occur in the corresponding eigenspaces. The 
nomenclature A\, A2, E, Fi, F 2 follows fll6f . 



jn/ J 


degeneracy 


irreducible representations of O 


-2 


5 


E®F 2 





3 


F 2 


2 


3 


F x 


4 


1 





The remaining 3-dimensional eigenspace of J with eigenvalue j max i = carries the 
3-dimensional self-representation of O and hence corresponds to the symmetric anti- 
ground state with 5 = which can be vizualized as a octahedron in spin space. Its 
energy is E max = 0. Again, as in the case of the cube, there is a weakly symmetric 
co-planar state with the same energy and theorem 1 yields E max (S) = ^fS 2 . This state 
has the form of a (not necessarily regular) hexagon such that opposing vertices of the 
octahedron in real space are mapped onto opposing vertices of the hexagon in spin space. 

The above RBS energy bounds can also be obtained more simply: Since Hq = 
S 2 — (Sf 6 + S 2 5 + S34), where = Si + Sj and J — 1, the minimal energy E min (S) is 
obtained for = 2 as E min (S) = S 2 — 12. Similarly the energy is maximal for = 
which yields E max (S) = S 2 (l - §) = |S 2 . 

6.4-3. The cuboctahedron The cuboctahedron is the quasi-regular polyhedron obtained 
by joining the mid-points of the cube's edges with their nearest neighbors, see 0, ||, 
and figure ||. 

It can be decomposed into 8 triangles with disjoint edges and it is 3-colorable, see 
figure 2. Hence there is a coplanar ground state with 5 = giving rise to an RBS- 
parabola 

E min (S) = J(^S 2 -2A), (89) 

since j — 4 and j m m = —2. The eigenvalues of J together with their degeneracy and 
irreducible representations of O are summarized in table 1. 

This has been calculated by using the formula connecting the characters with the 
multiplicity of irreducible representations, c. f. | 16f , 4.2.31b. From this table it is obvi- 



ous that the trivial irreducible representation A\ corresponds to the total anti-ground 
state, E is spanned by the complex eigenvector of the coplanar ground state, and the 
self-representation Fi corresponds to the cuboctahedron in spin space with the energy 
Emaxify — 12 j max i = 24 J. It remains to identify the symmetrical states which generate 
the two 3-dimensional irreducible representations F2 corresponding to the eigenvalues 
and —2 of J. 
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The vertex vectors r^fi = 1, . . . , 12 of the cub octahedron may be represented, 
up to normalization, by integer 3-vectors with components —1,0, 1 and exactly one 0- 
component. Then a state vector s M may be defined by the following rule: Invert the 
first component of after and set the other two components to 0. Here "after" is 
understood cyclically, e. g. (1,1,0) i— ► (—1,0,0). Thus we obtain a state, see figure [?| 
where the spins s /t of adjacent vertices are orthogonal. Indeed, this state corresponds 
to the 3-dimensional eigenstate of JJ with eigenvalue which is transformed under O 
according to F 2 . 

A second state s' is obtained by a permutation of the r M . Each is mapped onto 
that vertex where the above-defined state vector s M points to: 

s' ^ = 2s^ — r^, (not normalized) (90) 

Alternatively, s'^ is obtained from by the following rule: Invert the first component 
of r M after and leave the other unchanged. For example, (1, 1,0) h- > (—1, 1,0). Also 
this state generates a 3-dimensional subspace of the eigenspace of J with j min = —1,1 
which transforms under O according to F 2 . 

Geometrically, s' is a figure in spin space formed by four Stars of David, see figure 
Thus the cuboctahedron has a coplanar ground state with S = as well as a non- 
coplanar one with the same energy. 

6.5. y -symmetrical systems 

We will consider the two remaining Platonic solids as well as the quasi-regular icosi- 
dodecahedron and again calculate the decomposition of the eigenspaces of the adja- 
cency matrix J into 3^-hreducible subspaces for these three cases. As in the case of 
C-symmetric systems, the eigenvalue j max i always corresponds to the 3-dimensional 
self- representation of y called However, we have not found coplanar states realizing 
jmaxi and hence the question whether E max (S) will be a parabola remains open. 

The situation for ground states is different from the previous examples: For the 
icosahedron and the dodecahedron we find symmetric ground states with S — 
corresponding to the 3-dimensional irreducible representation F 2 of y, but no coplanar 
ones. Moreover, numerical calculations yield the results 

E m7n nar ~ -43.0614 . . . J > E mm = -20V5J for the dodecahedron, (91) 

E mln nar « - 24 J > E mm = -12V5J for the icosahedron. (92) 

Hence we conjecture that there are no coplanar ground states and hence, according to 
theorem 2, E m i n (S) will not be an exact parabola. Of course, even overwhelming nu- 
merical evidence cannot be considered as a rigorous proof. We can only strictly exclude 
symmetric coplanar ground states, since there are no two-dimensional irreducible repre- 
sentations of y. This is in accordance with proposition |] and the fact that the recently 
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Table 2. Eigenvalues j n of the adjacency matrix J of the icosahedron (1st column) 
together with their degeneracy (2nd column). In the 3rd column the irreducible 
representations of y are indicated which occur in the corresponding eigenspaces. The 
nomenclature A, Fi,F2,G,H follows Jit}] . 



jn /J 


degeneracy 


irreducible representations of y 


-VE 


3 


F 2 


-l 


5 


H 




3 


Fi 


5 


1 


A 



Table 3. Eigenvalues j n of the adjacency matrix J" of the dodecahedron (1st column) 
together with their degeneracy (2nd column). In the 3rd column the irreducible 
representations of y are indicated which occur in the corresponding eigenspaces. The 
nomenclature A, Fx, F 2 , G, H follows Q. 



jn/J 


degeneracy 


irreducible representations of y 


-VE 


3 


F 2 


-2 


4 


G 





4 


G 


1 


5 


H 


y/E 


3 


F 


3 


1 


A 



discovered 5 = ground states of the icosidodecahedron |§ are coplanar and weakly 
symmetric, but not symmetric. 

For the two ^-symmetrical Platonic solids we will investigate more closely the 
3-dimensional geometry of the ground states. By numerical simulation of a heat 
bath at zero temperature, Ch. Schroder has found the angles between adjacent 
spins to be ai ~ 116.6° for the icosahedron and an ~ 138.2° for the dodecahedron. 
From the knowledge of the corresponding irreducible subrepresentation of the natural 
representation of y in Mr we can now exactly determine the symmetric ground state, 
see proposition 0. It turns out that these states are well-known stellated geometrical 
structures, called "Great Icosahedron" for the ground state of the icosahedron and 
"Great Stellated Dodecahedron" for the ground state of the dodecahedron, see [0 and 
||. Thus the angles between adjacent spins are just the angles between neighboring 
vertices of these stellated structures and hence have the exact values 

c\r = arccos(— ^), ap = arccos(— ^), (93) 

in agreement with the numerical findings of [pOfl . 
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Table 4. Eigenvalues j n of the adjacency matrix J of the icosidodecahcdron (1st 
column) together with their degeneracy (2nd column). In the 3rd column the 
irreducible representations of y are indicated which occur in the corresponding 
eigenspaces. The nomenclature A, F\, F%, G, H follows Jl6| . 



jn/J 


degeneracy 


irreducible representations of y 


-2 


10 


H®H 


1 - V5 
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F 2 


-1 


4 


G 
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4 


G 


2 
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1 + V5 


3 


Fi 


4 


1 


A 



7. Summary 

We have tried to give a systematic survey of the properties of (relative) ground states of 
classical Heisenberg spin systems with particular emphasis on symmetric systems. Fur- 
ther we have devised various methods of ground state construction, e. g. extension of 
local ground states, construction of symmetric ground states from irreducible represen- 
tations of the system's symmetry group, and the construction of relative (anti-) ground 
states from coplanar ground states with S = 0. The latter procedure yields upper and 
lower parabolas as the boundaries of the S'-resolved energy spectrum. Thus we have a 
sufficient condition for a system to be parabolic. Moreover, under certain assumptions, 
this condition can also be proved to be necessary. 



The above issues are also relevant for the quantum theory of Heisenberg spin 
systems. It has been shown for various cases |IU|],|2"I|] that the shape of the S'-resolved 
energy spectrum of the quantum system is very well approximated by the curves E min (S) 
and E max (S) of the corresponding classical system, if only the individual spin quantum 
number s exceeds some moderate value, say s = 2. Hence one can predict semi- 
quantitative features of the quantum spectrum if one knows the classical spectrum and 
s>2. 

For some quantum systems, like the icosidodecahedron with s = 5/2, where the di- 
agonalization of Heisenberg Hamiltonian is totally impractical, approximation methods 
such as the density matrix renormalization group method (DMRG) are currently able 
to provide estimates for the lower boundary of the eigenvalue spectrum. In fact, for 
the icosidodecahedron the estimates |22| are in very good quantitative agreement with 
the lower boundary of the parabolic spectrum of the classical system. This agreement 
suggests that the DMRG results are also close to those for the quantum system with 
s = 5/2. For other classical systems with the same symmetry group y, viz. the dodec- 
ahedron and the icosahedron we expect non-parabolicity, although the corresponding 
quantum systems may exhibit approximate rotational bands. 
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In this sense our results and case studies are to a great extent relevant also for real, 
quantum spin systems. 
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Figure 2. Energy e versus square of total spin, S 2 , for the triangle with J\ — 0, J2 — 
1, J3 = —1 according to (|66|). The broken lines represent the bounding parabolas of 
theorem |[ 
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Figure 3. The shaded region represents the spectrum E versus S 2 of coplanar states 
of the pentagon as determined by numerical methods. The straight lines are the 
boundaries of the full spectrum. Also a part of the boundary of the coplanar spectrum 
which can be analytically calculated is displayed. 




Figure 4. Pentagonal star as an example of a spin system without coplanar ground 
states. 
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Figure 5. 3-dimensional ground state of the pentagonal star: The vertices correspond 
to 6 unit vectors in spin space and the edges correspond to the 10 bonds of the 
pentagonal star of figure EL see the remarks after (10). 
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Figure 6. The cuboctahedron (shaded figure) is obtained by joining the mid-points 
of the cube's edges with their nearest neighbors. 
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Figure 7. A stationary state of the cuboctahedron with energy E = 0. The spin 
vectors s M are shown as small arrows attached to vertices of the cuboctahedron. 
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Figure 8. A ground state of the cuboctahedron formed by 4 Stars of David. 



